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ABSTRACT 



We present the results of a new survey of 23 molecular clouds for the Zeeman effect in 
OH undertaken with the ATNF Parkes 64-m radio telescope and the NRAO Green Bank 
43-m radio telescope. The Zeeman effect was clearly detected in the cloud associated 
with the H II region RCW 38, with a field strength of 38 it 3 fiG, and possibly detected 
in a cloud associated with the H II region RCW 57, with a field strength of —203 it 24 
/xG. The remaining 21 measurements give formal upper limits to the magnetic field 
strength, with typical la sensitivities < 20 /xG. 



^Present Address: Deputy Vice-Chancellor, James Cook University, Townsville, QLD 4811, Australia 
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For 22 of the molecular clouds we are also able to determine the column density of 
the gas in which we have made a sensitive search for the Zeeman effect. We combine 
these results with previous Zeeman studies of 29 molecular clouds, most of which were 
compiled by Crutcher (1999), for a comparsion of theoretical models with the data. 

This comparison implies that if the clouds can be modeled as initially spherical 
with uniform magnetic fields and densities that evolve to their final equilibrium state 
assuming flux-freezing then the typical cloTid is magnetically supercritical, as was found 
by Crutcher (1999). If the clouds can be modeled as highly flattened sheets threaded 
by uniform perpendicular fields, then the typical cloud is approximately magnetically 
critical, in agreement with Shu et al. (1999), but only if the true values of the field 
for the non-detections are close to the 3a upper limits. If instead these values are 
significantly lower (for example, similar to the la limits), then the typical cloud is 
generally magnetically supercritical. 

When all observations of the Zeeman effect are considered, the single-dish detection 
rate of the OH Zeeman effect is relatively low. This result may be due to low mean field 
strengths, but a more realistic explanation may be significant field structure within the 
beam. As an example, for clouds associated with H II regions the molecular gas and 
magnetic field may be swept up into a thin shell, which results in a non-uniform field 
geometry and measurements of the beam-averaged field strength which are significantly 
lower than the true values. This effect makes it more difficult to distinguish magnetically 
subcritical and supercritical clouds. 

Subject headings: ISM: magnetic fields — ISM: molecules — ISM: clouds — ISM: 
kinematics and dynamics — ISM: individual (ROW 38, ROW 57) 

1. Introduction 

Magnetic fields are widely believed to influence the stability and dynamics of molecular clouds 
and to explain the supersonic line widths seen in all interstellar molecular spectral lines (see reviews 
by Heiles et al. 1993; McKee et al. 1993; and references therein). Yet these central ideas in molecular 
cloud and star formation physics are supported by remarkably few measurements of magnetic field 
strengths, due to the difficulty in measuring the Zeeman effect in molecular lines such as the 18 cm 
lines of OH. Excluding masers, for which the velocity dispersion and gas density are very uncertain, 
there are only 17 clouds in which the Zeeman effect has been detected above 3a (Crutcher 1999; 
Sarma et al. 2000; Crutcher & Troland 2000). 

As discussed in detail in the reviews by Heiles et al. (1993) and Crutcher (1994), the only 
viable technique available for measuring magnetic field strengths in molecular clouds is the Zeeman 
effect in spectral lines. The magnetic field reveals itself via the Zeeman effect as small frequency 
shifts, Avzj in the right and left circularly polarized (RCP and LCP, respectively) components of 



-3- 



the spectral line with respect to the frequency in the zero field case, Vq- The magnetic field strength 
can be determined by measuring this frequency shift. However, under most astrophysical conditions 
^ Azv, where Au is the full width at half maximum of the spectral line (the exceptions are 
masers, see e.g., Reid &: Silverstein 1990), and so detecting the shift between the RCP and LCP 
components due to the the Zeeman effect is difficult, and complete information about the magnetic 
field direction and magnitude is not obtainable. In this situation only the determination of the 
line-of-sight component of the field strength. Bigg (= BcosO), and its sign (i.e., toward or away 
from the observer) is possible. 

The most succussful efforts have been those of Crutcher and collaborators (e.g., Crutcher 
&; Kazes 1983; Kazes & Crutcher 1986; Troland, Crutcher & Kazes 1986; Crutcher, Kazes &: 
Troland 1987; Kazes et al. 1988, Goodman ct al. 1989; Crutcher et al. 1993; Roberts, Crutcher & 
Troland 1995; Crutcher et al. 1999a), mainly through observations of the OH transitions at 1665 
and 1667 MHz. Almost all successful detections of the Zeeman effect have used this technique, 
with observations of H I and CN also yielding success (e.g., Roberts et al. 1993; Brogan et al. 1999; 
Crutcher ct al. 1999b). The major telescopes that have been used are Nancay (OH), the Hat Creek 
85-ft (H I), the NRAO 43-m (OH, H I), Arecibo (OH, H I), the VLA (OH, H I), and the IRAM 
30-m (CN). 

All the reliable detections, and a number of upper limits were compiled by Crutcher (1999) 
for comparison with theory. In addition to the field strength or upper limit for each cloud in his 
analysis (27 clouds), Crutcher also compiled information on the column density, line width, number 
density, temperature and size for each molecular cloud. By comparing the data against a model 
spherical cloud with uniform density and uniform magnetic field, with corrections for the simple 
virial terms, Crutcher concluded that (i) internal motions are supersonic, (ii) the ratio of magnetic 
to thermal pressure is sufficient for magnetic fields to be important, (iii) the mass-to-magnetic flux 
ratio is about twice critical, which suggests that static magnetic fields alone may be insufficient 
to provide cloud support, (iv) the kinetic and magnetic energies are approximately equal, which 
implies that MHD waves and static fields are equally important in cloud dynamics, and (v) magnetic 
field strengths scale with gas densities in agreement with the predictions of ambipolar diffusion, 
\B\ oc /9°'^^. However, we note that the scatter in this last result is significant, and it does not take 
into account measured upper limits to Bios- 

Shu et al. (1999) compared the data of Crutcher (1999) with a model of a highly flattened 

molecular cloud (Allen & Shu 2000), and found that the mass-to-magnetic flux ratio is approxi- 
mately critical. Considering the uncertainties in the measurements and the simple nature of the 
model, Shu et al. concluded that it could not yet be decided whether the clouds are generally 
magnetically subcritical (slow evolution, dominated by ambipolar diffusion) or supercritical (fast 
evolution, dominated by collapse and fragmentation), though the distinction is clearly very impor- 
tant. 

Further, the role of the magnetic field in cloud support has been questioned since magneto- 
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hydrodynamic (MHD) simulations indicate that unforced MHD turbulence decays in about a free- 
fall time (MacLow ct al. 1998; Padoan & Nordland 1999; Stone, Ostriker & Gammie 1998), and 
since examination of stellar ages indicates that star-forming molecular clouds may not require 
long-lived support from either turbulence or magnetic fields (Elmegreen 2000, Hartmann 2001). 

Since progress in this field clearly requires more observations, we have undertaken OH Zeeman 
observations of molecular clouds in the relatively unexplored southern hemisphere and observations 
of a sample of northern hemisphere molecular clouds not previously studied. 

2. Observations 
2.1. Parkes 

OH Zeeman observations of southern hemisphere molecular clouds were undertaken in 1995 
July and 1996 October with the Australia Telescope National Facility (ATNF)^ Parkes 64-m ra- 
diotelescope, using a HEMT receiver equipped with dual linear polarization feeds. This system 
provided cold-sky system temperatures of Tgys ~20-25 K. The OH ground state transitions at 
1665.40184 (hereafter, the OH 1665 line) and 1667.35901 MHz (hereafter, the OH 1667 line) were 
observed simultaneously in one bandpass of 4 MHz centered on 1666 MHz with a channel separation 
of 488.28 Hz (0.087 km s'^). The outputs from the linear feeds were crossed in a hybrid coupler 
located immediately after the low noise amplifiers of the receiver to produce right- and left-hand 
circular polarized outputs (RCP and LCP respectively). These signals were then passed through 
a double-pole, double-throw polarization switch to minimize the effects of instrument polarization 
and gain differences on the spectra, switching every ~10 seconds. The center bandpass frequency 
was switched by ±0.25 MHz every 4 minutes. Both senses of circular polarization were observed 
and recorded simultaneously. The FWHM beam size at these frequencies is ^--^13' and the main 
beam efficiency is ^0.7. 

To determine the sense of polarization in the output from the correlator, OH masers having 
strong circular polarization were observed at the start of each run (Caswell et al. 1980; Caswell &; 
Haynes 1983, 1987a). We used Zeeman observations of the deep OH absorption line toward Orion B 
to test our setup. In both observing periods we were able to reproduce the results of Crutcher & 
Kazes (1983; hereafter CK83) after ~5 hr integration time. By combining the 1665 and 1667 OH 
results from 1995 and 1996 a field strength of 35 it 2 //G is inferred for the subcomponent identified 
by CK83 as the source of the Zeeman pattern, essentially independent of hour angle. The Stokes / 
and V spectra for Orion B observed in 1996 are shown in Figure 1. These values are in agreement 
with the value determined by CK83 (38±1 /iG) using the Nangay telescope. It should be noted that 
the magnetic field strength inferred for Orion B is not unique, but depends on the properties (line 



^The Australia Telescope is funded by the Commonwealth of Australia for operation as a National Facility managed 
by CSIRO 
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strength, line width) of the subcomponent, which are not fully constrained by the observations. 

The displacement of the right- and left-circularly polarized beams on the sky, commonly called 
"beam-squint" (e.g., Troland &; Heiles 1982) can result in a spurious Zeeman effect if there exists 
a sufficiently large component of velocity gradient in the spectral line along the same direction 
as the squint. This effect is diminished considerably with the combination of an altitude- azimuth 
telescope such as Parkes and a long on-source integration time, since the source then rotates within 
the beam and the true Zeeman signal is modulated. For Parkes we measured a squint of ~15" for 
both observing runs, but with different position angles, most likely caused by the replacement of 
the focus cabin at the end of 1995. While this is a relatively large value for the squint, the fact that 
we were able to reproduce the results for Orion B both times with the same derived magnetic field 
strength suggests that the squint is not a problem for our observations. If a spurious Zeeman effect 
had been caused by squint then the inferred magnetic field strength for the OH 1665 line would be 
1.7 times greater than the OH 1667 line, which we do not see in our results. 



2.2. Green Bank 

OH Zeeman observations of northern hemisphere molecular clouds were undertaken in 1996 
September and 1998 June with the National Radio Astronomy Observatory (NRAO)^ Green Bank 
43-m telescope, using a dual channel cooled HEMT receiver with circular polarization feeds. Cold 
sky system temperatures of 20-25 K were recorded with this system. The OH 1665 and 1667 MHz 
lines were observed simultaneously in both polarizations with frequency switching with a channel 
separation of 1220 Hz (~0.22 km s"^). To minimize the effect of instrument polarization and gain 
differences on the spectra polarization switching was performed at a rate of ~1 Hz. The setup is 
similar to the one used by Crutcher et al. (1993). The FWHM beam size at these frequencies is 
~18' and the main beam efficiency is ~0.7. 

As was the case for the Parkes observations, we tested the setup with OH Zeeman observations 
of Orion B. The Zeeman effect was clearly seen during both observing sessions, and by combining 
the 1665 and 1667 OH results for 1996 and 1998 we infer a magnetic field strengths of 34 ± 3 fiG 
for the subcomponent identified by CK83. 

We measTircd the beam squint of the 43-m telescope to be no more than 1". Beam squint 
is more important for an equatorially mounted telescope such as the 43-m compared to an alt-az 
telescope, if the squint happens to align with a velocity gradient within the source being observed. 
However, the chances of such an alignment axe small, and the change in hour angle will not cause 
any change in alignment between the squint and velocity gradient. 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 
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2.3. Mopra Observations 

As OH and ^^CO are believed to trace similar column densities (e.g., Myers et al. 1978; Heiles 
et al. 1993), we have mapped a number of the southern clouds (those observed at Parkes) in 
the ^^CO J = 1 — ^ transition, in order to provide an independent measurement of the column 
density, and to resolve more structure than is possible with the poor angular resolution provided 
by the OH observations. The observations were undertaken in April 1997 with the 22-m ATNF 
Mopra radiotelescopc. At the ^'^CO J = 1 ^ frequency of 110201.37 MHz only the inner 15-m 
of the surface was illuminated at the time of our observations, providing an effective beam size 
of 45". The receiver was a dual linear polarization SIS receiver with a total system temperature 
measured to be ~150 K in good weather. The 16384 channel autocorrelator was configured for two 
IFs, providing a spacing between channels of 0.17 km s~^ over a velocity range of ~220 km s~^ 
(1024 spectral channels with a 64 MHz bandwidth) at the ^^CO J = 1 ^ frequency. One IF 
was centered on the ^^CO J=1^0 frequency, the second being centered on the SiO 1^ = 1, 
J = 2 ^ 1 transition at 86243.44 MHz for regular pointing checks using SiO masers listed in the 
SEST Handbook (http://www.ls.eso.org/lasilla/Telescopes/SEST/SEST.html). 

The mapping was performed on a 1' grid in position switching mode, with an on-source integra- 
tion time of 4 minutes per position. Removal of first-order baselines was sufficient in most cases to 
correct for any system effects not removed by the reference observation. The spectra were calibrated 
offline for atmospheric attenuation and the variation of telescope beam efficiency with elevation, 
and the intensity scaled to the scale of SEST by observations of Orion (aigso = 05'^32"^47!0, 
Ji950 = -05°24'23") in the ^"'CO J = 1 transition, assuming r^(Orion) = 13 K (Appendix of 
the SEST Handbook). The conversion to Tmb is given by T^b = T\/rj^^'^ where r]^^'^ = 0.7 near 
110 GHz. In general we attempted to map down to at least the 50% level of the peak emission in 
each source. 



3. Source Selection 

3.1. Southern sources 

Our initial search for suitable sources centered on the OH studies of Robinson, Goss & Manch- 
ester (1970), Goss, Manchester & Robinson (1970), Manchester, Robinson & Goss (1970), Robinson, 
Caswell & Goss (1971), Caswell & Robinson (1974), and Turner (1979). The spectral resolution 
of these previous observations, many undertaken with Parkes, was generally insufficient to detect 
weak, narrow maser features, and so we reobserved all our potential targets with the setup de- 
scribed here. Wc also reobserved a number of other sources from these lists with a greater velocity 
coverage, as well as some sources from Dickel & Wall (1974), Toriseva, Hoglund &: Mattila (1985), 
and Chan, Henning & Schreyer (1996). None of the sources we surveyed from these lists showed 
any significant OH absorption that was not previously known. 
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In order to select suitable Zeeman candidates we required observations with substantial im- 
provements in sensitivity and spectral resolution over the earlier OH surveys of southern HII regions. 
We thus rc-observed ~80 sources drawn from these lists to assess their suitability. To determine our 
primary candidates we used the "sensitivity estimator" given by Troland (1990; see also Goodman 
1989) to estimate the integration time required per source to reach our target sensitivity of 3a ~30 
fiG. In order to have the greatest sensitivity to the Zeeman effect, the spectral line profiles of the 
selected sources should be narrow and strong against a weak continuum source (so that the system 
temperature is as low as possible), and free of maser emission. Using these criteria we chose the 
9 sources for our primary candidate list, selecting those sources with the lowest integration times 
required to reach the target sensitivity. 

Previously unreported masers discovered in the survey are discussed in Appendix B. 

3.2. Northern sources 

Crutcher and co-workers have previously observed the most promising OH Zeeman candidates 
in the northern sky (Crutcher 1999). However, a number of promising sources remain to be studied. 
The best list to date is the extensive survey of OH by Turner (1979) using the NRAO 43-m telescope. 
Two of the previous detections of the Zeeman effect, S88B (Crutcher et al. 1987) and S106 (Kazes 
et al. 1988), have only modest OH line absorption depths of ~ 0.5 as observed with the 43-m, 
but they are among the few sources with clearly detected Zeeman patterns. The Turner catalog 
contains a number of sources with similar line absorption depths. To select sources for Zeeman 
observations with the 43-m, we re-observed a number of candidates from the Turner catalog, as 
well as sources identified as H II regions from the IRAS database, sensitive continuum surveys, and 
water maser surveys (e.g., Chan et al. 1996; Carpenter, Snell & Schoerb 1990; Wouterloot, Brand 
& Fiegle 1993; Wouterloot & Brand 1989). We observed ~200 sources, and identified 11 Zeeman 
candidates using the same criteria as for the Parkes sample. 

The final source list is given in Table 1. Column 1 lists the source name, columns 2-3 the 

source position, and columns 4-6 the results of fitting a Gaussian to the OH 1667 line profile. 
Column 7 lists the continuum source temperature, where appropriate, and column 8 list the total 
on-source integration time for the Zeeman observations. 

4. Results & Analysis 

In Zeeman experiments, the magnetic field reveals itself as frequency shifts in the right and left 
circular polarized components {Ircp and Ilcp^ respectively) of the spectral line. A full treatment 
of the Zeeman effect under astrophysical conditions may be found in Sault et al. (1990; see also 
Heiles et al. 1993), and we include only a brief discussion of the analysis here. In the situation where 
the frequency shift is much smaller than the line width (which is always the case for non-maser 
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lines), the difference or Stokes V spectrum, V = Ircp — Ilcp, is only sensitive to the hne-of-sight 
component of the magnetic field and can be modeled as 

where h is the "Zeeman factor" and equals 3.27 and 1.96 Hz /xG~^ respectively for the 1665 and 
1667 MHz lines of OH, B is the magnetic field strength, is the angle between the field lines and the 
line-of-sight, (3 is the "gain term" arising from gain differences between the RCP and LCP signal 
paths (Sault et al. 1990) and I is the Stokes / spectrum {=Ircp + Ilcp)- The Zeeman term in 
equation (1), {b/2)Bcos9{dI /du), reveals itself in the V spectrum as the characteristic sideways "S", 
or "Zeeman pattern". The line-of-sight field strength Bios = BcosO is inferred from equation (1) by 
least-squares fitting the right hand side of this equation to the observed V spectrum. Other fitting 
techniques are discussed and analysed in Sault et al. (1990). 



4.1. Fitting dl/di^ to V 

In Table 2 we present the results of our Zeeman analysis. In this Table column 1 lists the source 
name, column 2 the velocity of the OH feature, column 3 the magnetic field strength determined 
from the OH 1665 transition with its one sigma uncertainty, column 4 as for column 3 but for 
the OH 1667 transition, and column 5 the weighted mean line-of-sight field strength Bio^ (found 
by combining the results for the 1665 &: 1667 lines) and its weighted one sigma uncertainty as 
obtained by standard error analysis. 

In the majority of sources the line profiles arc simple and well fitted by a single Gaussian 
profile. Notable exceptions are ROW 38, ROW 57, W7, and GlO.2-0.3. In these cases, alterna- 
tive approaches to performing the Zeeman analysis are described in CK83, or Heiles (1988) and 
Goodman Sz Heiles (1994). In the method of CK83, the V spectrum shows a clear "S" structure 
signifying the Zeeman effect, but it is too narrow to arise from the entire / line profile. By inte- 
grating the V spectrum and fitting the result with a Gaussian it is possible to infer the properties 
{Visr, AV) of the component of the / spectrum responsible for the "S". The reality of such a com- 
ponent is strengthened if the same properties are derived for both the OH 1665 and 1667 lines. 
CK83 applied this method to Orion B (see Fig. 1), and we have successfully used this technique 
to analyse ROW 38 (sec below). In Heiles (1988) and Goodman &; Heiles (1994) the / profile is 
broken down into independent Gaussian subcomponents, and the V spectrum is then modeled as 
the sum of derivatives of these components, with unique values of Bios assigned to each component. 
We have applied this method to those sources for which a single Gaussian could not adequately be 
fitted to the / spectrum and where the CK83 method is inappropriate, because there is no clear 
"S" in the V spectra. These are indicated in the notes to Table 2. For GlO.2-0.3, we were unable 
to fit a consistent set of Gaussians to the OH 1665 and 1667 lines. 

In order for the derived magnetic field strength Bios to be considered as a detection we require 
that the field strengths measured in both OH lines be consistent with each other, | -Biges — -Bi667 | 
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< (Tfi, and the field strength be greater than SaB, \ Bi^s \ > Sas- Only 2 sources, ROW 38 and 
RCW 57, can be considered as detections using these criteria. These are discussed below. We 
consider the remaining 21 measurements to give upper limits to Bios averaged over our 13' (Parkes) 
and 18' (Green Bank) beams (see also §5.1). 

Details of other properties for the individual sources, such as column densities, are given in 
Appendix A (see also Table 3). 



4.2. Determination of Column Densities 

The intensity of a spectral line from a uniform medium in front of a background region may 
be expressed as 

T^b = f[J{Tea,)-JiTcm-eM-r)] , (2) 

where T^b is the main-beam brightness temperature, r is the optical depth of the cloud in the molec- 
ular transition, / is the beam filling factor, Tex is the excitation temperature of the transition, Tc 
is the area-weighted mean background temperature (the cosmic background temperature T^g plus 
any contribution from a background continuum source, T5), and J{T) = {hiy/k)/(exp(hi'/kT) — 1) 
is the Planck function. For sources observed in absorption against a compact continuum source, / 
represents the fraction of the beam area covered by the continuum source. At the frequencies of the 
OH ground state transitions near 1.6 GHz J(T) = T for the temperatures observed in molecular 
clouds (T < 100 K), and so 

T,,b = f[Tex-Tc][l-eM-r)] ■ (3) 

The observed antenna temperature is proportional to T^fe, with a proportionality constant -qmb 
(the main beam efficiency). With this knowledge it is possible to determine the column density of 
OH, iV(OH), which is given by (Goss 1968; Turner & Heiles 1971) 

where v is the frequency of the transition and A is its Einstein A coefhcient, gi is the statistical 
weight for level i and Qu is the statistical weight for the upper state of the transition. Thus, for the 
OH 1667 line 

Ar(OH) = 2.25 X 10- (Zh^) ,,,,, em- . (5) 

For clouds unassociated with H II regions, previous observations have shown that T^x is gen- 
erally about 6 K (e.g.. Turner 1973; Crutcher 1977, 1979), and we assume this value here. Clouds 
associated with H II regions are warmer, and so we assume Tex = 10 K for those clouds. For a 
cloud seen in absorption against a typical H II region Tq S> Tex and Tex drops out of equation (3) . 
We measure the product fTc at the telescope as T^j^s (on-source) - T^j^^ (off-source), and the op- 
tical depth is then determined from equation (3), and the column density from equation (5). To 
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determine the molecular hydrogen column density we assume a standard OH abundance ratio of 
[H]/[OH] = 2.5 X 10'' (Crutcher 1979), which is uncertain by about a factor two. High-angular res- 
olution observations indicate that the OH abundance in unshockcd gas found in PDRs associated 
with H II regions could be enhanced by as much as a factor of 5 over this value (Roberts et al. 
1995; Sarma et al. 2000). Comparison of the data for NGC 6334 in Table 3 suggests that this is 
not a concern for the single-dish observations. Sarma et al. (2000) use an OH abundance five times 
greater than the dark cloud value and determine iV(H2) ~ 5 x 10^^ cm^^, whereas we find iV(H2) 
~ 2 X 10^^ cm~^. The discrepancy would be greater if we used the same abundance as Sarma et al. 
(2000). Further support for the OH abundance ratio we use is given below. 

Where it is possible to estimate the ^^CO column density independently, or it has been previ- 
ously determined (see appendix for individual cases), the values of the molecular hydrogen column 
density derived independently from OH and ^^CO generally agree to within a factor two. For 
sources where a map of the ^^CO emission exist, the line profile of ^^CO averaged over the map is 
remarkable similar to OH, with similar line velocities (Fig. 5). The agreement in line profile and 
the molecular hydrogen column density derived from ^^CO and OH observations supports the as- 
sumption that these molecules trace similar physical conditions and the values for Tg^ and [H]/[OH] 
used here. 

4.3. RCW 38 

The Stokes / and V spectra for RCW 38 are shown in Figure 2. The characteristic "S" shape 
of the Zeeman effect is clearly evident in the V spectra, but is not well fitted by a scaled version of 
dl/du. As discussed above, this suggests that the observed Zeeman profile is not due to the entire 
/ profile. The Zeeman profile is due to a subcomponent of the observed line, and we integrated the 
V spectra to determine its properties. The integrated V spectra are well fitted by single Gaussians, 
essentially with identical properties for the OH 1665 and 1667 lines. The line velocity and width 
for this component is given Table 2. From the fit of the derivative of this Gaussian to the OH lines 
we infer a value for the magnetic field of Bi^s = 38 it 3 fiG. 

RCW 38 was observed during both observing sessions at Parkes. The results are consistent 
with each other. This helps to validate our result as a true Zeeman detection, as the direction of 
beam squint for the two observing sessions was different, even though the magnitude of the squint 
was large for Zeeman observations. In addition, no velocity gradients are apparent in a map of 
the OH absorption made during the 1996 observations. Corroborating evidence for the reality of 
the RCW 38 Zeeman detection comes from high angular resolution interferometric observations 
of the OH absorption made with the ATNF Compact Array (Bourke et al. 2001 in prep.). These 
observations indicate a compact, spatially isolated absorption feature with line velocity and width 
similar to those of the sub-component we infer as giving rise to the Zeeman effect. 

RCW 38 is a bright H II region which has recently been observed in the near-infrared with 
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the VLT (Alves et al. 2001 in prep). In a separate paper we present a detailed molecular line 
and continuum study of the region (Bourke et al. 2001 in prep.) which indicates that the Zeeman 
effect arises from a molecular clump on the western edge of the H II region, possibly associated 
with the bright 10 //m peak IRSl (Smith et al. 1999). Assuming T^.^ = 10 K we determine rieer 
= 0.15, iV(OH) = 7.6 X IQi^ cm'^, and iV(H2) = 9.3 x lO^^ cm-^. ^^CO J =1^0 observations 
imply Ar(H2) = 4 X 10^^ cm~^ (assuming [H2]/[^^C0] = 7 x 10^), lower than derived from the OH 
observations, which may suggest that the OH is sampling denser gas (e.g., Roberts ct al. 1995 for 
S106). Observations of CS J = 2 ^ 1 and C^^O J = 1 ^ support this view (Bourke et al. 2001 
in prep.). 

4.4. RCW 57 

The Stokes / and V spectra for RCW 57 are shown in Figure 3. The / profile is well fitted 
by two overlapping Gaussians, a deep absorption component at Visr = -26.1 km s~^ and a weaker 
component at Visr = -22.6 km s~^. The ^^CO J = 1 ^ channel maps shown in Figure 4 reveal 
two spatially distinct clouds. The cloud in the north-east corresponds in velocity to the blue-shifted 
OH component, and the larger cloud to the south-west to the red-shifted OH component, which is 
the component for which we claim the detection of the Zeeman effect. The sum of ^^CO emission 
spectra over the mapped region is remarkably similar to the OH absorption spectra (Fig. 5), and 
is well fitted by two Gaussians with similar line velocities and widths as those fitted to the OH 
spectra. 

The formal results for RCW 57 indicate a solid detection of the Zeeman effect in the —22.6 
km s""^ component, with Bios = —203 it 24 /xG. However, inspection of the V spectra in Figure 3 
raises some doubt about this result. The blending of the two components makes it difficult to say 
with certainty that the Zeeman effect has been detected. At velocities greater than about —22.2 
km s~^ there is no overlap between the OH velocity components, and the Zeeman effect is clearly 
seen, but only part of the "S" is evident, due to the line overlap at lower velocities. As the ^^CO 
observations show two spatially distinct clouds, follow-up OH Zeeman observations with the ATNF 
Compact Array should be able to resolve this question. 

Located at a distance of 3.6 kpc (kinematic, e.g., Caswell &: Haynes 1987b) or 2.4 kpc (Persi 
et al. 1994, based on the possible association of HD 97499 with the molecular cloud), RCW 57 
(NGC 3576) is an H II region with ongoing star formation (e.g., Persi et al. 1994). The few published 
studies of this region suggest it is a "typical" H II region with associated molecular cloTids. Whiteoak 
(fe Gardner (1974), observed the 5 GHz transition of H2CO in absorption toward RCW 57 with 
a 4?2 beam, and found two velocity components with almost identical properties to the Gaussian 
components we fitted to our OH data, suggesting that at least part of the OH absorption arises in 
somewhat denser gas (10^ cm~^) than is usually sampled by OH (10^ cm~^). For the component 
at -22.6 km s"^ we find nmr = 0.06 and iV(OH) = 8.9 x 10^^ cm-^ or iV(H2) = 1.1 x lO^^ cm-^ 
compared with iV(H2) = 1.3 x 10^^ cm~^ from ^^CO Mopra observations. This is the component 
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where we have a probable detection of the Zeeman effect, with Bi^g = —203 ± 24 /xG. For the 
component at -26.1 km we find neey = 0.2 and iV(OH) = 1.5 x 10^^ cm~^, or A^(H2) = 
1.9 X 10^^ cm~^ compared with N{B.2) = 1-5 x 10^^ cm~^ from ^^CO Mopra observations. 



5. Discussion 

Of the 23 Gaussian components we could identify along the line-of-sight toward the 20 sources 
observed, the Zeeman effect was clearly detected in one source (RCW 38) and possibly in one 
component of another (RCW 57), as indicated in Table 2. In Table 2 we label our 23 components 
as those directly associated with H II regions (those clouds whose OH absorption line velocities are 
similar to the recombination line velocity of the background H II region) and those associated with 
(cold) molecular clouds along the line-of-sight. We find that 2 of 13 components associated with 
H II regions are detected, while none of the 10 components associated with cold molecular clouds 
are detected. The median la sensitivities for the H II and cold-cloud samples are respectively 11 
and 9 /xG. 

A comparison of OH optical depth against line width shows no clear trend, nor does a com- 
parison of integrated OH line area against optical depth. The molecular clouds showing Zeeman 
detections do not have properties clearly different to the other clouds in our sample. 



5.1. The magnetic field and virial equilibrium 

If the gas is well coupled to the magnetic field, then the large scale magnetic field can provide 
support against gravitational collapse. An important parameter in the discussion of support from 

large scale magnetic fields is the magnetic flux-to-mass ratio, ^/M. To examine whether this large 
scale field alone can provide support to molecular clouds, wc ignore terms in the virial theorem 
other than the magnetic (Ai) and gravitational (W) energies, and assume that these energies are 
in equilibrium. This allows us to determine a critical mass, Mcr, when gravitational forces are 
balanced by magnetic stresses (this has been discussed in detail in McKee et al. (1993); see also 
McKee 1999). It is found that M = W implies 

where G is the gravitational constant, and c$ is a numerical factor which depends on the geometry 
under consideration. For a spherical "parent" cloud with uniform density and uniform magnetic 
field, c$ ~ 0.12 (Mouschovias & Spitzer 1976; Tomisaka, Ikeuchi & Nakamura 1988) for the final 
equilibrium configuration, which is flattened along the field lines (assuming flux- freezing) . If the 
cloud is an isothermal disk (i.e., sheet-like) with a constant flux-to- mass ratio (uniform field and 
uniform column density), then c$ = l/27r 0.16 (Nakano & Nakamura 1978; Shu et al. 1999). 
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The critical mass can be related to observable quantities with measurements of the magnetic 
field strength B and column density N (McKee et al. 1993; Crutcher 1999), so that 

$ B 

— = (7) 

where m is the mean particle mass (2.33mH for 10% He). The column density in equation (7) 
is the total column density, A'^ = 1.2A^(H2), for 10% He. The magnetic field strength is the 
total field strength, whereas Zeeman experiments where the Zeeman shift is less than a line width 
are only sensitive to the field along the line-of-sight, and hence measure Bi^g. As discussed in 
Crutcher (1999), we can statistically estimate the line-of-sight field strength for a large ensemble of 
measurements of uniform fields oriented randomly with-respect-to the line-of-sight, with the result 
that < Bios > = I 1/2 and < Bios^ > = \ B |^/3. In addition, for the sheet-like geometry a 
correction factor to the observed column density is required, again due to the orientation with- 
respect-to the line-of-sight. In particular, if the magnetic field is preferentially perpendicular to the 
sheet, then the correction factors to the observed column density Ngijs and Bi^g are not independent. 
Following the same argument as for the ensemble of magnetic field measurements given above, we 
find that B/N = {Bios/Nohs){cos^e) for one line-of-sight, so that B/N = < Wios/Nobs > for the 
average over a distribution of equally likely line-of-sight directions. In terms of the critical magnetic 
flux- to-mass ratio, 
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where the subscript "n" indicates that the values of <I>/M have been normalized by the critical 
value (<I>/M)cr, and F is the geometric correction factor, which is 2 for the initially uniform sphere, 
and 3 for the isothermal sheet. Evaluating the constants leads to 
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= 4 X 10^^—^ cm^ iiG (sheet) . (10) 
N{E.2) 



The values of (<I>/M)n for the different geometries considered above are given in Table 3. In 
this Table we also include the data for B and N compiled by Crutcher (1999) and recent Zeeman 
results from Sarma et al. (2000; NGC 6334) and Crutcher & Troland (2000; L1544). Therefore we 
believe Table 3 gives the most complete summary available of Zeeman measurements in extended 
molecular clouds (excluding small-scale OH maser emission). In column 1 is the source name, 
columns 2 & 3 give the line-of-sight field strength Bios and its log, and column 4 lists the log of 
the molecular hydrogen column density. Column 5 lists the values of ($/M)n assuming an initially 
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uniform sphere. Column 6 lists the values of ($/M)n assuming a sheet-like geometry. The values 
in Table 3 are shown graphically in Figure 6. In this figure are plotted the measured field strength 
B 

los against the molecular hydrogen column density A'^(H2), and overlayed with loci of ($/iVf)n for 
the uniform sphere model (a) and the sheet-like model (b). 

Considering the values listed in Table 3 and plotted in Figure 6, the results for the spherical 
parent cloud model are consistent with the clouds being slightly magnetically supercritical (($/M)n 
< 1), the same result as found by Crutcher (1999), who only considered the detections he tabulated. 
The addition of our large number of upper limit measurements to his sample strongly reinforces this 
result. As noted by Crutcher (1999), assuming an initially uniform spherical cloud with a uniform 
magnetic field may be an oversimplification. In particular, if clouds are supported primarily by 
static magnetic fields with uniform direction they will become flattened as matter flows inward 
along field lines, and the field will become pinched along the direction perpendicular to the initial 
field direction. Therefore, the results for the sheet-like cloud may be more realistic. Considering 
only those clouds where the Zeeman effect was detected. Table 3 shows that the results for this 
geometry are consistent with ($/M)n ~ 1, and the clouds are in approximate equilibrium between 
magnetic pressure and gravity, the same result found by Shu et al. (1999), though they did not 
distinguish between the detections and the upper limits. Because of the statistical corrections 
applied to B and N, and the uncertainties in measuring A^, a more definitive conclusion cannot be 
reached at present. 

If we assume for the clouds where the Zeeman effect was not detected that the true magnetic 
field strengths are equal to the 3a upper limit values (listed in Table 3) , then for the sheet- like model 
the mean flux-to-mass ratio is <($/M)n> = 1.3 ±0.9 1 (excluding G20.8-0.1 and G29.9-I-0.0 as 
the sensitivity of these observations is significantly less than the other data), the same results as 
for the detections. If the magnetic field within the beam is uniform, the true values are not likely 
to be greater than 3a for any of the clouds. On Figure 6 we also indicate the la limits for the 
non-detections. Assuming a uniform field, if the true field strengths are closer to the la limits than 
to the 3(7 limits then for the sheet-like model the mean flux-to-mass ratio is <{^/M)n> = 0.4±0.3, 
significantly lower than the critical value, which implies that in general the clouds are magnetically 
supercritical by a factor 2-3. As discussed below, field structure within the beam could reduce the 
flux-to-mass ratio by a factor 3-4 over its true value. However, beam dilution will also result in 
measurements of the column density lower than the true value for the gas in which the magnetic 
held is studied, which will increase the fiux-to-mass ratio. 

It is important to note that a number of the sources located at the upper end of Figure 6 (log 
A^(H2) > 22.6) are due to high angular resolution VLA Zeeman observations. In almost all of these 
cases the initial detection of the Zeeman effect was made with low-angular resolution single-dish 
observations. Consequently, weaker field strengths by a factor ~3-4 were reported compared to the 
subsequent higher-angular resolution observations which reveal the field structure. Therefore there 
is an inconsistancy between the lower and upper ends of Figure 6, as the lower end contains many 
upper limit values from low-angular resolution observations, and the upper end contains many 
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clouds where the Zeeman effect has been detected in higher-angular resolution inteferometric and 
single-dish observations. Observations at higher angular resolution of sources showing only upper 
limits to B in low- angular resolution single-dish observations are required. If the lack of detection 
of the Zeeman effect is due mainly to field structure within the large single-dish beam, then these 
higher resolution observations should in many cases reveal the Zeeman effect (e.g., Brogan et al. 
1999 for M17; Sarma et al. 2000 for NGC 6334) and allow us to explore more fully the high column 
density range of the plot and the structure in the field. 

Table 3 also shows that there is no clear example of a magnetically subcritical cloud, as has 
been noted by Nakano (1998). The cloud toward RCW 57 where we have a possible detection of 
the Zeeman effect, and hence a measurement of B, has a high value of ($/M)n for the geometries 
considered here, implying that it is in a magnetically subcritical state. As discussed earlier, high 
resolution Zeeman observations of RCW 57 are required to examine the validity of this result. 

5.1.1. Non-uniform Fields Associated with H II Regions 

The foregoing results show that the flux-to-mass ratio deduced from the observations may or 
may not be magnetically critical, depending on whether the magnetic field is modeled as threading 
a region whose geometry is spherical or planar. Here we point out another geometrical effect which 
can influence the conclusion of magnetic criticality. The field direction may not be uniform within 
the telescope beam, particularly for a situation common in single-dish Zeeman observations: the 
molecular gas is associated with an H II region, and is seen in absorption against the H II region. 
In such a case, the gas and field around the protostar may be compressed into a thin shell by the 
pressure of the H II region, so that the field is largely tangent to the surface of the H II region. The 
similarity of magnetic and H II region pressures has been noted for the well-studied source W3 OH 
(Reid, Myers &: Bieging 1987). The resulting geometry differs from the cases discussed above, since 
the field direction is not uniform, but lies tangent to a layer surrounding the H II region, except 
near the "poles" which mark the original field direction. A sketch of this situation is shown in 
Figure 7. 

If the shell is perfectly spherical, a uniform field B with initial angle (p from the line-of-sight can 
be shown to have a line-of-sight component, averaged over the forward hemisphere, of < Bios > = 
{2/3)B cos (f); and the mean over a random distribution of initial angles is then ^ Bios ^ = B/3. 
Thus the true field strength is a factor of 3 greater than the typical line-of-sight component. Also, 
this geometrical factor should increase further if weighting due to the varying brightness across 
the face of the H II region is taken into account. Then the brightest part of the H II region, at its 
projected center, is where the field direction is most nearly perpendicular to the line-of-sight, and 
the faintest part of the H II region, at its projected edge, is where the field lies most nearly along 
the line-of-sight. 

If this geometrical picture is correct, one should expect to see evidence for the corresponding 
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field structure in higher-resolution Zeeman maps and in submiUimeter polarization maps of H II 
regions. The "polarization holes" noted in many submiUimeter maps of H II regions (Hildebrand 
et al. 1995) merit study as possible counterparts of the magnetic polar regions described here. 

In the shell model the total column density of the gas in the forward hemisphere is unchanged 
from before to after the formation of the shell, so the flux-to-mass ratio requires the same correction, 
by a factor ~3, as does B alone. However caution is required in applying the results of the shell 
model to obtain an estimation of the flux-to-mass ratio, since the present-day field is in the plane 
of the shell, rather than perpendicular to it; and since the original flux-to-mass ratio may have 
changed during the expansion of the H II region, if the expansion induced some flow along fleld 
lines. Furthermore, as indicated by the beam-filling factors in Table 1, in low-angular resolution 
single-dish observations the source is generally much smaller than the beam. This may significantly 
reduce the effect of the geometry discussed here in accounting for the low detection rate. 

6. Conclusions 

We have presented the results of a survey of 23 molecular clouds for the Zeeman effect in OH. 
For 22 of these clouds we were also able to determine the column density. We combined our data 
with the data for 29 clouds analysed by Crutcher (1999), and compared the combined results with 
simple cloud models. 

The data in Table 3 and Figure 6 are generally consistent with a flux-to-mass ratio less than 
its critical value by a factor of a few, if the model cloud is initially a uniform sphere threaded 

by a field of uniform strength and direction, independent of whether the true field strengths for 
the non-detections are close to or significantly less than the 3cr upper limits. Such a magnetically 
supercritical cloud should be highly unstable to collapse and fragmentation. These conclusions are 
in agreement with those of Crutcher (1999). 

If instead the model cloud is a highly fiattened sheet threaded by a uniform perpendicular 
field, the data are generally consistent with a critical flux-to-mass ratio, implying a critically stable 
system (assuming that the true values of B for the non-detections are close to the 3a upper limits) . 
These conclusions are in agreement with those of Crutcher (1999) and those of Shu et al. (1999). 

However, if the true values of B for the non-detections are closer to the la values, then for the 
sheet-like model cloud the data are generally consistent with a flux-to-mass ratio significantly less 
than critical, implying that the typical cloud is significantly magnetically supercritical. 

When both negative and positive detections of the Zeeman effect are considered, the single- 
dish detection rate of the OH Zeeman effect is relatively low, less than 10% in the present survey. 
This low rate may be due simply to low mean field strengths, which is the simple assumption we 
have used in our analysis. A more realistic explanation of the low detection rate may be a selection 
effect which tends to decrease the field strength inferred from a single-dish Zeeman absorption 
observation of the most common survey target, an H II region. Expansion of an H II region may 
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compress the gas and field into a shell-like configuration, decreasing the Zeeman eflFect compared 
to that from a uniform field of the same magnitude and mean direction. Interaction with the H II 
region may also lead to line-of-sight reversals of the field (e.g., Brogan et al. 1999), which tend to 
reduce the Zeeman effect by a factor of a few when observed with coarse resolution. 

These considerations imply that the flux-to-mass ratio of the typical cloud associated with an 
H II region may appear critical or somewhat supercritical on the size scale to which single-dish 
observations are sensitive. But it is difficult to infer this ratio accurately from such low-resolution 
observations and from models of uniform fields, because the field is likely to have important unre- 
solved structure. 

In summary the principal results and conclusions of this study are: 

1. The Zeeman effect was detected in a molecular cloud associated with the H II region 
RCW 38, with a field strength of 38 ± 3 //G. 

2. The Zeeman effect was possibly detected in a cloud associated with the H II region RCW 57, 
with a field strength of -203 ± 24 /xG. 

3. If the molecular clouds arc modeled as initially uniform density spheres with uniform 
magnetic fields, then generally they are magnetically supercritical, and therefore unstable to grav- 
itational collapse. 

4. If the molecular clouds are modeled as highly flattened isothermal sheets, then those with 
detections of the Zeeman effect are approximately magnetically critical. Those clouds with non- 
detections of the Zeeman effect (and hence only upper limit values to the magnetic field strength) are 

also approximately magnetically critical, but only if the true field strengths are close to the 3a upper 
limits reported here. If instead the true field strengths are signifcantly lower, for example equal to 
the la upper limits, then on the average the clouds are approximately magnetically supercritical 
by a factor 2-3. 

5. For clouds associated with H II regions, the molecular gas may be swept up into a thin 
shell, resulting in a non-uniform magnetic field geometry. When this geometry is observed at low 
angular resolution, as is the case with single-dish observations, field cancellation within the beam 
will occur, resulting in measurements of the field strength which are significantly lower than the 
true values. 

6. A number of upper limits to the field strength exist for clouds observed in single-dish 
experiments. Those observations only sample moderate column densities and the clouds should be 
observed at higher angular resolution to examine whether the lack of a detectable Zeeman effect is 
due to changes in field direction within the single-dish beam. Such observations may hold the most 
promise for further studies of magnetic field strength in local molecular clouds. 
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A. Individual Clouds 
A.l. Parkes 

RCW 38 - see § 4.3 

Carina Molecular Cloud - Associated with t] Carina. It has been mapped in OH by Dickel 
&; Wall (1974), who found two peaks in the absorption, closely corresponding to the dark lanes 
either side of the optical emission associated with r) Carina. The position we have observed is the 
western peak. Brooks et al. (2001 in prep) are currently undertaking a multi-wavelength study of 
the Carina Molecular Cloud, and preliminary results from this study have been reported (Brooks, 
Whitcoak &: Storey 1998). The distance has been determined by Tapia et al. (1988) to be 2.5 kpc. 
Assuming T^x = lOK we determine rise? = 0.04 and A''(OH) = 5.2 x lO^'* cm~^, which imply N(il2) 
= 6.6 X 10^^ cm~^, the same result we obtain from the CO and ^^CO data of Brooks et al. (2001 
in prep). Dickel & Wall (1974) find nmr = 0.033 and A^(OH)/re^ = 4.9 x 10^^ cm'^ K'^ at the 
same position. 

Chamaeleon I- Chamaeleon I is a well studied molecular cloud (see e.g., the review by Schwartz 
1991) forming mainly low-mass stars, located at a distance of 160 pc (Whittet et al. 1997). Though 

it is easier to detect the Zeeman effect in strong absorption lines from molecular clouds associated 
with H II regions, there is a clear need for more measurements of the magnetic field strength 
in cold molecular clouds such as Chamaeleon (Crutcher et al. 1993), which is why it has been 
included in this study. There is evidence for a large scale ordered magnetic field in the region 
(McGregor et al. 1994), and this together with the bright OH lines (at least for a thermal emission 
source) makes Chamaeleon I a prime Zeeman candidate among dark clouds. The cloud has been 
mapped in OH with the Parkes radiotelescope by Toriseva et al. (1985) and the extinction has been 
determined throughout the entire cloud by Cambresy et al. (1997) using near-infrared colours. We 
assume = G K and so nee? = 0.48, iV(OH) = 6.3 x 10^^ cm'^, and iV(H2) = 7.8 x 10^^ cm'^. 
For comparison, the NANTEN and DENIS data in the same region (Hayakawa et al. 1999) imply 
iV(i3C0) = 5 - 10 X 10^5 cm-2, or iV(H2) = 3.5 - 7 x lO^^ cm-^. At the position we observed, 
Toriseva et al. (1985) obtain Ar(OH) 7 x 10^^ cm'^. 



-19- 



RCW 57 -see § 4.4 

G326.7+0.6 - Very little is known about this region. The radio continuum source is probably 
associated with the optical H II region RCW 95. The OH component at -44.8 km s~^ has the same 
velocity as the H109a recombination line (Caswell & Haynes 1987b), and is associated with the H II 
region at a distance of 3.0 kpc (kinematic). The narrow line width of the -21.6 km s^^ component 
suggests that it arises from an unassociated dark cloud along the line-of-sight to the continuum 
source, with a kinematic distance of 1.5 kpc. For the -44.8 and -21.6 km s~^ components we 
determine Ar(OH) = 9.5 x 10^^ cm~^ and 3.1 x 10^^ cm~^ respectively, assuming Tg^ = 10 and 6 
K, respectively. These values imply N(}i2) = 1.2 x 10^^ cm~^ (compared to A^(H2) = 2.3 x 10^^ 
cm-2 from ^^CO observations) and A^(H2) = 3.8 x lO^^ cm'^ (Ar(H2) = 2 x lO^^ cm'^ from ^^CO) 
respectively. 

G327.3-0.5 - Possibly associated with RCW 97 at a kinematic distance of ~3.3 kpc, the only 
detailed study of the molecular gas associated with G327.3-0.5 has been reported by Brand et al. 
(1984) who mapped the cloud in the ^^CO J = 2 ^ 1 transition, and found a massive ('-^lO^ -^o) 
cloud associated with a number of infrared sources (Frogel Sz. Persson 1974), and a possible bipolar 
molecular outflow associated with the southern-most group of infrared sources. Circularly polarized 
masers are present in both OH lines, and unfortunately extend well into the blue side of the main 
OH absorption line. For the absorption line we find A^(OH) = 2.3 x 10^^ cm~^, assuming T^x = 
10 K, which implies A^(H2) = 2.8 x 10^^ cm~^, comparing favourably with that determined from 
the Mopra ^^CO observations of iV(H2) = 3 x 10^^ cm . The Zccman effect is clearly seen in 
three of the masers present in the 1667 line. However, since the properties of the molecular gas 
associated with masers is poorly understood we simply state that the magnetic field strengths we 
infer from the observed splitting (Reid & Silverstein 1990) are +1.3, +2.9, and +7 mC, for masers 
with velocities at —80, —68 and —52 km s~^, respectively, though the identification of the Zeeman 
pair for the —52 km s~^ component is not certain. 

G343.4-0.0 - Another southern H II region which has not been studied, G343. 4-0.0 is not 
clearly associated with any RCW region. The two main absorption features are at —27 km s~^ 
(near the recombination line velocity) and the narrow, deep line near 6 km s~^. Other weak 
absorption lines are seen, but the signal-to-noise (S/N) is < 50 so they are not considered in the 
Zeeman analysis below. For the same reason the —27 km s~^ feature is not used in the Zeeman 
analysis. We somewhat arbitrarily assign the near kinematic distance of 2.8 kpc to the H II region. 
Unlike the sources discussed above, G343. 4-0.0 has not even been included in many of the major 
molecular line and maser surveys of the southern galactic plane. Whiteoak & Gardner (1974) 
observed the 6 km component in formaldehyde absorption, and OH absorption observations 
were reported in the extensive OH survey by Turner (1979). The 6 km s~^ cloud is clearly local, 
and may be associated with the Lupus clouds, which have a similar velocity and lie about 10° above 
the plane at a similar galactic longitude. It may also be associated with the 6 km s^^ component 
seen toward NGC 6334 and NGC 6357 (W22; Kazcs & Crutchcr 1986; Crutcher et al. 1987; Massi, 
Brand & Felli 1997) near longitude 351°. Kazes & Crutcher (1986) detected the Zeeman effect in 
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the 6 km component toward NGC 6357 (W22B) with a magnetic field strength of — 18=b 1 ^G, 
while Crutcher et al. (1987) claim a detection at a nearby position (W22A) with a field strength 
of —32 lb 9 ;uG. The large scale CO maps made with the Columbia 1.2 m radiotelescopc (e.g., 
Bronfman et al. 1989; Bitran et al. 1997), suggest an association between all the clouds near 6 
km s~^. If associated with Lupus, a distance <200 pc is appropriate. For this velocity component 
we find N{0}i) = 4.3 x 10^^ cm~^, which implies N{B.2) = 5.4 x 10^^ cm"^ compared with N{B.2) 
= 2 X 10^^ cm~^ derived from the ^^CO observations. 

NGC 6334 - NGC 6334 is a well studied H II region located at 1.7 kpc (Neckel 1978). Visually 
it is very prominent, with numerous nebulous spots, and forms a fine pairing with the nearby NGC 
6357 (W22). Important studies include the CO observations of Dickel, Dickel & Wilson (1977), the 
OH observations by Brooks (1995; published in Brooks &: Whiteoak 2001) and the recent PhD thesis 
by Kraemer (1998), who studied [O I], [C II], CO, CS and ammonia throughout the star-forming 
molecular "ridge". There are two main absorption components, one at —4 km associated with 
the H II region, and one at 6 km s~^. As discussed above for G343.4-0.0, the 6 km s~^ component is 
local and covers a large area, but has not been studied in any detail and its distance and properties 
are not well known. Recently Sarma et al. (2000) presented Zccman observations of the —4 km s^^ 
component in OH and H I made with the VLA. Following the nomenclature of Rodriguez, Canto 
&; Moran (1982), they measured the Zeeman effect in OH against source A {Bios= 152 ± 17 //G) 
and in H I against source D (-93 ± 13 //G), source E (-180 ± 29 //G), and NW of D (169 ± 55 
/xG). The position we observed covers source D and the position NW of D. Sources A and E lie 
just outside our primary beam. The ATCA OH channel maps presented by Brooks (1995) show 
strong absorption against source D, and only weak absorption against source E and NW of D, 
which suggests that the absorption we observe in the single-dish observations is mainly due to 
source D. We do not detect the Zeeman effect in OH. It is clear that the magnetic field changes 
direction within the Parkes beam, as indicated by the signs of Bios for D and NW of D, so the lack 
of detection may be due to field reversal and therefore cancellation within the beam. From our 
observations we find A/'(OH) = 1.4 x 10^^ cm~^, which implies A^(H2) = 1.8 x 10^^ cm~^. From 
our limited ^^CO observations we find A^(H2) = 1 x 10^^ cm~^, while Kraemer &: Jackson (1999) 
obtain a value of A^(H2) = 2.2 x 10^^ cm^^ toward source D from their CO observations. 

Troland and co-workers (Troland 1995, private communication) have searched for the Zeeman 
effect toward NGC 6334 with the Nancay telescope and report a 3a upper limit for the 6 km s^^ 
component of ~20 /LtG. For this cloud we find iV(OH) = 1.2 x 10^^ cm'^ and iV(H2) = 2.8 x lO^i 
cm~^, assuming Tgx = 6 K. Prom ^^CO observations we obtain A^(H2) = 2 x 10^^ cm~^. 

G8. 1+0.2 - The recombination line velocity for the G8.H-0.2 H II region is ~20 km s~^, which 
implies a kinematic distance of ~3.5 kpc. The spectra show two velocity components near 13 and 17 
km s"^ respectively, and a circularly polarized maser pair near 22 km s~^. The S/N ratio of the 13 
km s~^ component is too low to be considered in the Zeeman analysis. For the 17 km component 
A^(OH) = 2.9 X 10^^ cm"2 ^^^^ iV(H2) = 3.6 x 10^^ cm'^. We have no ^^CO observations for an 
independent check of these values. The Zeeman effect is clearly seen in the maser lines, from which 
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we can infer a magnetic field strength in the masing region of ~ +3 mG. No detailed study of this 
H II region has been reported. 

A.2. Green Bank 

W7 ~ Also known as 3C 123, the absorption lines are due to the Taurus Molecular Cloud 
(Crutchcr 1977). Bregman et al. (1983) attempted to detect the Zeeman effect in H I using the 
WSRT, and reported a 3a upper limit of 16 /nG. From our observations we determine rise? = 0.09 
and Ar(OH) = 2.4 x 10^^ cm~^. In the same region Wouterloot & Habing (1985; their region f), 
find iV(OH) ^8x lO^^ cm-^. 

NGC 2264-IRS2 - Also known as IRAS 06382+0939, it was observed during periods when the 
other sources on our list were not available. The region has been observed in C^^O J = 1 — >^ by 
Schwartz et al. (1985), who determine iV(H2) f« 1 x 10^^ cm-^^ and in ^^CO J = 2 ^ 1 by Wilking 
et al. (1989), who find A'"(H2) = 2 x 10^^ cm~^. Prom our OH observations we obtain A^(H2) = 

5.3 X 10^^ cm~^. IRS2 is associated with the northern molecular core, while IRSl (Allen's infrared 
source) is found in the southern core. Due to the weak OH emission lines our limit to Bios is not 
as low as has been determined for other dark molecular clouds (e.g., Crutcher et al. 1993). 

GlO.2-0.3 - Associated with the W31 region at a distance of ~2.0 kpc. As the Zeeman effect 
was not detected and we were not able to decompose the complex line profile with the same set of 
Gaussians for both OH transitions, no further analysis was performed. 

G14.0-0.6 - The OH line velocity is similar to that of the recombination line (Lockman 1989), 
with a kinematic distance of ~2.4 kpc. It may be associated with the other sources observed here 
with velocities near 20 km s-^ including IRAS 18153-1651 and G14.5-0.6/IRAS 18164-1645 (Jaffe, 
Stier & Fazio 1982). No studies of the molecular gas in this region arc known. We find A?"(H2) = 

1.4 X 10^^ cm~^. The continuum source may be associated with IRAS 18151-1707. 

IRAS 18153-1651 - Located in the same region as G14.0 0.6, a similar distance is assumed. 
Jaffe et al. (1982) detected it as a far-infrared source (G14.21-0.53) between 40-250 nm, and 
from their ^^CO observations we infer a column density of A'^(H2) = 5.3 x 10^^ cm~^. From our 
observations we find A^(H2) = 1.3 x 10^^ cm~^. 

G14.5-0.6/IRAS 18164-1645 - Our Green Bank observations of these two sources overlap. 
They are part of the same group of sources located around 2.4 kpc (which includes G14.0-0.6 and 
IRAS 18153-1651) and were observed by Jaffe et al. (1982) in ^^CO (G14.5-0.6 = G14.33-0.64 and 
IRAS 18164-1645 = G14.43-0.69). Using their observations of ^^CO we infer A^(H2) = 6.4 x 10^^ 
cm-2 for G14.5 0.6 and iV(H2) = 5.7 x 10^2 cm'^ for IRAS 18164-1645. From our OH observations 
we obtain iV(H2) = 1.7 x 10^2 cm-^ for G14.5-0.6 and A^(H2) = 1.2 x 10^^ cm'^ for 18164-1645. 
We know of no other studies of these clouds. 

G20. 8-0.1 - The H II region is associated with IRAS 18264-1652, and has a recombination line 
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velocity of ~56 km s~ . The absorption line we see in OH is unassociated with the H II region, 

with a velocity of 6.8 km s~^. This velocity implies a kinematic distance of only 600 pc, which is 
unreliable, but docs indicate that the molecular cloud is probably nearby. No studies of this cloud 
were found in the literature. From our OH observations we infer A^(H2) = 1.2 x 10^^ cm~^. 

G29. 9+0.0 - The background continuum source is the very well studied UCH II region G29.96- 
0.02 (e.g., Pratap, Megeath &; Bergin 1999), with a recombination line velocity near 100 km s~^. 
As is the case for G20. 8-0.1, our OH observations arise from a foreground region near 8 km s~^, 
in an anonymous molecular cloud. Its position and velocity suggest that it is associated with the 
Aquila Rift series of clouds (Dame et al. 1987), which includes the well know Serpens star-forming 
region at a distance of ~300 pc (De Lara, Chavarn'a-K, & Lopez Molina 1991). This may also 
be true for G20.8-0.1. Prom our OH observations we find A^(H2) = 2.1 x lO^^ cm'^ along this 
line-of-sight. 

G78.5+1.2 - Also known as L889, OH Zeeman observations were reported by Crutchcr ct al. 
(1993), who found a 3a upper limit of 6 //G, much lower than our limit of 30 /xG. We observed a 
position ~17' away from their position. At this position we find A'^(H2) = 3.9 x lO^-*^ cm~^, while 
Crutcher et al. (1993) find N{'H.2) = 10^^ cm~^ at the position they observed. Dickel, Seacord 
&; Gottesman (1977) observed 2 and 6 cm H2CO absorption toward this region and derived a 
maximum column density of N(H.2) = 6 x 10^^ cm~^. 

G80.9-0.2 - Associated with the compact source DR22, which has a recombination line velocity 
near km s~^, implying a kinematic distance of 3 kpc. The OH lines are near 7 km s"-'', for which 
no molecular line observations have been published. Prom the OH observations we determine A'^(H2) 
= 7.8 X 10^1 cm-2. 



B. Masers 

G312.1+0.3 - (R.A., Dec.) = 14:05:06, -60:56:09 (B1950.0). No maser emission is evident in 
the OH 1665 line. In the OH 1667 line a single component is seen in both RCP and LCP within 
the absorption line. This maser has Visr = —46.1 km s~^ and AV = 0.42 km s~^. 

G327.3-0.5 - see Appendix A.l. The feature at —80 km s~^ has not previously been reported. 

G331.4-0.0 - (R.A., Dec.) = 16:07:14.0, -51:23:26) (B1950.0). A maser feature at Visr = 
—62 km with two velocity components is present in the OH 1665 line in both RCP and LCP 
(absorption is present at Visr = —46 km s~^). In the OH 1667 line weak (T| ~ 0.5 K) maser 
emission may be present in the absorption line in both RCP and LCP, with Visr = —44.6 km s"^ 
and AV = 0.43 km s^^ 
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(^^Not well represented by a Gaussian 



(b) Emission line 
(=)IRAS 06382+0939 
(■^^ Beams overlaps 

(^)L889, nearby position observed by Crutcher et al. 1993 

Note. — For sources where the electron temperature (Te) has been measured from radio recombination lines 
observations, the source filling factor / can be estimated. The average is ~6500 K, which implies a range in 
/ of 0.002-0.013. Similar results are found for sources where a size estimate is available. 
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Table 2. Magnetic Field Strengths 



IN dililc 


T/; 
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5 ±18 


26± 12 


y 
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-43 ± 23 


-25± 17 


y 
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-16 ±25 
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6.4 
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'^'^^ Method of CK83 used to infer properties of Gaussian subcomponent from integrated 
V spectrum. The line width of this Gaussian is AV = 2.3 km s~^ 

('^^Method of Heiles (1988) and Goodman & Heiles (1994) used, as the Gaussian com- 
ponents overlap. The Zeeman fitting is performed simultaneously on both velocity com- 
ponents. 

('^^The 1665 MHz result is listed in parentheses, as interference prevents a meaningful 
fit. 

(^)Two velocity components. Only one component is fitted, as the other is too narrow. 

(•^^ Complex line not well represented by Gaussian components. Results reported here 
are for the whole line profile. 
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Tablo 3. MolocTilar Cloud Analysis 
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21.73 


< 0.3 


< 0.7 




NGC 6334 (vl) 


< 12 


< 1.08 


22.25 


< 0.1 


< 0.3 


(c) 


NGC 6334 (v2) 


< 12 


< 1.08 


21.46 


< 0.8 


< 1.6 


(c) 


G8.1+0.2 


< 36 


< 1.56 


22.56 


< 0.2 


< 0.4 




1 J.T'OOn 1 < * A T 1 \r 

VJlCCll JDclliK 














W7 


< 18 


< 1.26 


21.50 


< 1 1 

\, J- • J- 


< 2.3 




N2264-IRS2 


< 48 


< 1.68 


21.73 


< 1.8 


< 3.5 




G14.0-0.6 


< 56 


< 1.75 


22.15 


< 0.8 


< 1.6 




IRAS 18153-1651 


< 51 


< 1.71 


22.10 


< 0.8 


< 1.6 




G14.5-0.6 


< 54 


< 1.73 


22.22 


< 0.7 


< 1.3 




IRAS 18164-1645 


< 51 


< 1.71 


22.08 


< 0.9 


< 1.7 




G20.8-0.1 


< 69 


< 1.84 


21.09 


< 11.2 


< 22.1 




G29.9+0.0 


< 78 


< 1.89 


21.33 


< 7.0 


< 14.0 




G78.5+1.2 


< 30 


< 1.48 


21.59 


< 1.5 


< 3.0 




G80.9-0.2 


< 27 


< 1.43 


21.89 


< 0.7 


< 1.4 




Crutcher 1999 














W3 OH 


3100 


3.49 


23.7 


1.2 


2.4 




DR 21 OHl 


710 


2.85 


23.6 


0.4 


0.7 




Sgr B2 


480 


2.68 


23.4 


0.4 


0.8 




M17 SW 


450 


2.65 


23.1 


0.7 


1.4 




W3 (main) 


400 


2.60 


23.2 


0.5 


1.0 




S106 


400 


2.60 


22.8 


1.3 


2.5 




DR 21 0H2 


360 


2.56 


23.3 


0.4 


0.7 
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Table 3 — Continued 



Name 


Bios 


Log B 


Log N 


($/M)n 


($/M)n 


Notes 






(mG) 


(iV(H2) cm-i) 


(Sphere) 


(Sheet) 




OMC-1 


360 


2.56 


23.2 


0.5 


0.9 




NGC 2024 


87 


1.9 


22.9 


0.2 


0.4 




S88 B 


69 


1.84 


22.3 


0.7 


1.4 




Bl 


27 


1.43 


21.9 


0.7 


1.3 




W49 B 


21 


1.32 


21.6 


1.1 


2.1 




W22 


18 


1.26 


22.2 


0.2 


0.5 




W40 


14 


1.15 


22. 


0.3 


0.6 




p Oph 1 


10 


1. 


21.7 


0.4 


0.8 




L1544 


11 


1.04 


22.3 


0.1 


0.2 


(d) 


NGC 6334 


150 


2.18 


22.7 


0.6 


1.2 


(e) 


0MC-N4 


< 300 


< 2.48 


23.1 


< 0.5 


< 0.9 


(f) 


Tau G 


< 7 


< 0.85 


21.6 


< 0.4 


< 0.7 


(f) 


T 1 fiQ 


< io 


< i.io 




< i.y 


^ 'X 1 

<. o.i 


lij 


L1647 


< 11 


< 1.04 


22.1 


< 0.2 


< 0.3 


(f) 


p Oph 2 


< 13 


< 1.11 


21.6 


< 0.7 


< 1.3 


(f) 


TMC-1 


< 12 


< 1.08 


21.9 


< 0.3 


< 0.6 


(f) 


L1495W 


< 9 


< 0.95 


21.6 


< 0.5 


< 0.9 


(f) 


L134 


< 10 


< 1. 


21.3 


< 1.0 


< 2.0 


(f) 


TMC-IC 


< 9 


< 0.95 


21.9 


< 0.2 


< 0.5 


(f) 


L1521 


< 9 


< 0.95 


21.7 


< 0.4 


< 0.7 


(f) 


L889 


< 6 


< 0.78 


22. 


< 0.1 


< 0.2 


(f) 


Tau 16 


< 7 


< 0.85 


21.7 


< 0.3 


< 0.6 


(f) 



(^)(vl) = -26.1 km component, (v2) = -22.6 km s~^. 
^^^(vl) = -44.8 km s^"*^ component, (v2) = -21.6 km s~"^. 
*^^)(vl) = -3.8 km component, (v2) = 6.3 km s~^. 
('i)Prom Crutcher & Troland(2000) 

(^)Prom Sarma et al. (2000). Corresponds to our NGC 6334 (vl). 

(^)Note that the values Hsted here for Bi^g are 3a b for consistency with our 
other results, rather than Bigg + 3(7 as in Crutcher (1999). 
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Fig. 1. — The Zeeman effect toward Orion B, as observed with the Parkes radio telescope in October 
1996. On the left is the data for the OH 1665 line, on the right the OH 1667 line. In (a) is the OH 
line profile (average of right- and left-circular polarizations; solid line) with the Gaussian component 
responsible for the Zeeman effect indicated (dotted line). In (b) is the Stokes V spectrum (right- 
minus left-circular polarization; light grey line), with the derivative of the Gaussian fitted to derive 
the magnetic field strength (continuous black line). The value of Bios inferred from the fits are 
indicated. 

Fig. 2. — The Zeeman effect toward RCW 38, as observed in October 1996. On the left is the 
data for the OH 1665 line, on the right the OH 1667 line. In (a) is the OH line profile (average of 
right- and left-circular polarizations; solid line) with the Gaussian component responsible for the 
Zeeman effect indicated (dotted line). In (b) is the Stokes V spectrum (right- minus left-circular 
polarization; light grey line), with the derivative of the Gaussian fitted to derive the magnetic field 
strength (continuous black line). The value of Bi^s inferred from the fits are indicated. The OH 
1665 and 1667 Gaussians are assumed to have the same line depth, and so are optically thick. 

Fig. 3. — The Zeeman effect toward RCW 57. On the left is the data for the OH 1665 line, on the 
right the OH 1667 line. In (a) is the OH line profile (average of right- and left-circular polarizations; 
solid line) with the Gaussian components whose sum best fits the line profile indicated (dotted and 
dashed lines). In (b) is the Stokes V spectrum (right- minus left-circular polarization). Maser 
emission is present in the OH 1665 V spectrum at velocities near —25 km s~^. The feature near 
— 27 km s^^ in the OH 1667 spectrum is probably due to gain differences between the right- and 
left-circular polarizations, though very weak maser emission cannot be ruled out. In (c) is the 
V spectrum (light spectrum) with the best overall fit (dark line), consisting of the sum of the 
derivatives of the Gaussian components shown in (a), and a gain term, which is a scaled version of 
the OH line profile. The maser emission in the OH 1665 V spectrum has been excluded when fitting 
the spectrum. The feature near —27 km s~^ in the OH 1667 spectrum appears to be well fitted 
by the gain term component in the overall fit. In (d) is shown the V spectrum (light spectrum) 
overlayed with the derivative of the Gaussian component at —22.6 km s~^, which is one component 
of the fits shown in (c). The value of the Bigg inferred from this component is indicated at the 
bottom of the panel. 

Fig. 4.— 13 CO J = 1 ^ channel maps for RCW 57. The central velocity of each panel is indicated 
in the upper right. The Mopra beam size is shown in the lower left corner of the upper left panel 
as a filled circle, and the Parkes OH beam size is indicated by the dashed circle, which is centered 
on the position observed for the Zeeman effect (indicated with a star). The greyscale ranges from 
1 to 12 K km s-i in TJ. The contour levels are 30, 50, 70, and 90% of 12 K km 8"^ 

Fig. 5. — Average of RCW 57 ^^CO J = 1 ^ spectra compared with the OH 1667 MHz absorption 
spectrum. The ^^CO spectrum is offset by 0.5 K and the OH spectrum by —0.5 K. The map of 
i^CO J = 1 ^ emission is shown in Figure 4, with the area covered by the OH observations. 
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Fig. 6. — Magnetic field strength Biog plotted against column density A'^(H2) for the sources ob- 
served in this paper (squares) and for previously published results (circles; Crutchcr 1999; Sarma 
et al. 2000; Crutchcr & Troland 2000). The large symbols represent detections of the Zccman effect, 
while the small symbols represent 3a upper limits. The lines drawn down from the upper limits 
represent the shift of the 3a upper limits to the la limits for each source. In (a) are shown the 
loci for different values of the magnetic flux-to-mass ratio, assuming an initial configuration of a 
uniform density sphere, as described in Crutchcr (1999). The values have been normalized to the 
critical value, so that values > 1 are subcritical, and < 1 are supercritical. In (b) are shown loci 
assuming an infinite sheet geometry ( "highly flattened molecular cloud" ) , like that discussed in Shu 
et al. (1999). 

Fig. 7. — Sketch of magnetic field configuration near a spherical H II region. The field lines are 
initially uniform, and retain their original distribution far from the H II region. As the H II region 
develops, it expands and compresses the surrounding molecular gas and magnetic field into a thin 
spherical shell. This configuration implies a low detection rate for Zeeman absorption observations 
made with a beam larger than the H II region. For any line-of-sight direction the line-of-sight 
field component, averaged over the forward hemisphere, is substantially less than the line-of-sight 
component of a uniform field having the same mean magnitude and direction. 
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